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Abstract 

We present a calculation of the loop-induced processes gg — > HH and gg —* HHH, and 
investigate the observability of multi-Higgs boson production at the CERN Large Hadron Col- 
lider (LHC) in the Standard Model (SM) and beyond. While the SM cross sections are too 
small to allow observation at the LHC, we demonstrate that physics beyond the SM can lead 
to amplified, observable cross sections. Furthermore, the applicability of the heavy top quark 
approximation in two- and three-Higgs boson production is investigated. We conclude that 
multi-Higgs boson production at the SuperLHC is an interesting probe of Higgs sectors beyond 
the SM and warrants further study. 



1 Introduction 



The understanding of electroweak symmetry breaking is a major goal of the physics programme at 
the upcoming Large Hadron Collider (LHC) at CERN. In the Standard Model (SM), electroweak 
symmetry breaking is mediated by the Higgs mechanism, which predicts a fundamental scalar par- 
ticle, the Higgs boson [1]. At the LHC, the discovery of the Higgs boson is feasible in the entire 
theory-compatible mass region [2,3]. Theoretically, a heavy Higgs boson with mass of order 1 TeV 
in the vicinity of the non-perturbative regime is still feasible [4,5]. Electroweak precision measure- 
ments, however, favour a relatively light SM Higgs boson not too far above the direct LEP exclusion 
limit of 114.4 GeV [6]. 

In order to establish the Higgs mechanism and confirm the Higgs sector of the SM, the discovery 
of a Higgs- like boson is not sufficient. In addition, the predicted couplings of fermions and gauge 
bosons to the Higgs boson, as well as the trilinear and quartic Higgs self couplings need to be con- 
firmed experimentally. While Higgs-fcrmion and Higgs-gauge boson couplings are measurable with 
accuracies of 10-40% at the LHC [7] and in many channels considerably better at the International 
Linear Collider (ILC) [8], the measurement of the trilinear and quartic Higgs self couplings, which 
are probed in double and triple Higgs boson production, respectively, are more challenging. For low 
Higgs masses not too much above the LEP limit (m H < 140 GeV), the largest rates are obtained 
with the decay channel H — > bb dominating in this mass regime. Unfortunately, at the LHC this 
search channel is not viable, because of an overwhelming QCD background. The ILC, however, 
would allow the measurement of the trilinear Higgs self coupling to a precision of 20-30% in the low 
Higgs mass regime [8-11]. We note that in this regime a photon collider promises an even better 
determination of the trilinear Higgs coupling [12-14]. For higher Higgs masses, the ILC production 
cross section decreases due to the reduced phase space. However, in this mass region, the vector 
boson pair decay channels open up and allow for leptonic signatures that can be separated from the 
backgrounds at the LHC [8, 15-19]. 

The dominant production mechanism for Higgs boson pairs at the LHC is gluon fusion. We 
note that Higgs production in hadronic collisions can also proceed through bottom quark fusion, 
bb — > nH , but in the SM the corresponding LO as well as NLO [20] cross sections are negligible. In 
the MSSM, however, enhanced Yukawa couplings can lead to comparable cross sections for gluon 
and bottom quark fusion [21]. 

The gluon fusion loop amplitude was first presented in [22] . For neutral Higgs boson pairs in the 
Minimal Supersymmetric Standard Model (MSSM) the top/bottom loop contribution was evaluated 
in [23,24]. Charged Higgs boson pairs including squark effects were studied in [25]. Although the 
SM cross section for triple Higgs boson production at the LHC can be expected to be small [9, 10], 
this expectation has to be verified through explicit calculation. Only recently, a full calculation of 
the process gg — > HHH appeared in the literature [26] and confirmed that SM cross sections are 
indeed too small to be observable at the LHC. However, it has to be stressed that multi-Higgs boson 
production rates are very sensitive to physics beyond the Standard Model (BSM) and should be 
scrutinised carefully in high energy experiments, as they might point to new physics at high energy 
scales that arc not directly accessible at the given collider [27,28]. A non-standard heavy quark 
that receives its mass via the Higgs mechanism does not decouple [29], and therefore leads to a 
non-vanishing contribution in heavy quark loop-induced processes. Furthermore, contributions of 
higher dimensional operators might alter the SM cross section considerably [30]. Also, in certain 
Little Higgs Models the Higgs pair production cross section is significantly different from the SM 
cross section [31]. 

In this paper, we present our calculation of double and triple Higgs boson production via gluon 
fusion. This provides an independent check of the recent calculation in [26] , which employed different 
computational methods and tools. To go beyond the findings of [26], we study higher dimensional 
operator effects on production rates relative to the SM. The corresponding couplings are in principle 
only restricted by unitarity constraints [32] . We also analyse amplification effects in supersymmetric 
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(SUSY) two-Higgs-doublet models (2HDMs). 

The paper is organised as follows. After briefly reviewing Higgs boson properties relevant for 
our investigation in Section 2, we describe our calculation in Section 3. In Section 4, SM and BSM 
results are presented and discussed. Section 5 gives a brief summary. 

2 Higgs boson properties 

The SM Higgs mechanism provides fundamental mass terms for massive vector bosons and fermions. 
The coupling strength of the Higgs boson is proportional to the mass (squared) of the interacting 
fermion (gauge boson). The SM Higgs boson self interactions, induced by the scalar potential 

v -$(.'.-£)'. (i) 

are also proportional to the Higgs mass squared. In unitary gauge one has 

V=^ + |^ + ^ (2) 

with 

A 4 = h v = —g- . 3 

Relation ([3]) between the trilinear and quartic Higgs self couplings is a genuine SM prediction. To 
establish the SM Higgs mechanism, it has to be verified experimentally. As already pointed out, the 
discovery of the Higgs boson and the measurement of its couplings to fermions and gauge bosons 
alone are not sufficient. 

In general, the Higgs self couplings change in extensions of the SM. By allowing for higher 
dimensional operators of the type 

fe=l v 7 

the constraint Q is relaxed. Magnitude and sign of A3 and A4 can be arbitrary up to constraints 
imposed by unitarity. In order to guarantee the stability of the vacuum, only the sign of the highest 
power of the Higgs field has to be positive. We note that the addition of singlet Higgs fields preserves 
relation Q , but may lead to invisible Higgs boson decays and diluted Higgs signals [33] . 

The ability to measure the Higgs self couplings depends on the size of multi-Higgs boson cross 
sections. As will be discussed below, SM rates are very small at the LHC (see also [26]). It is 
thus interesting to consider extensions of the SM that allow for amplified event rates. For Higgs 
pair production in gluon fusion this has been studied in [23] in the context of the MSSM. Two 
amplification sources have been identified. Firstly, the top and bottom Yukawa couplings are altered 
due to the mixing of the Higgs fields. In the MSSM one has, at tree level, 
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Here, v± (i> 2 ) is the vacuum expectation value of the Higgs doublet with weak hypercharge — ^ (+^)- 
For sufficiently large tan/3, the bottom-loop contribution to the cross section becomes sizable in 
comparison to the top-loop contribution, leading to a larger production rate. Secondly, internal Higgs 
propagators can become resonant thereby enhancing the production rate. Both effects generally play 
a role in 2HDMs, as will be discussed and quantified for triple Higgs boson production below. 

We will now briefly review the features of 2HDMs that are important for our purposes. The 
general potential of the 2HDM is given by [34, 35] 

V($i,$ 2 ) = m^^^i + m^^^- (mi 2 $l$2 + h.c.) + 
+ y(<i> t 1 <i>i) 2 + Y(^ 2 ) 2 + 
+A 3 ($ t 1 $i)(^$ 2 ) + A 4 ($ t 1 $ 2 )(^*i) + 

+ {y( < J >t i < f > 2) 2 + [A 6 ($ t 1 <fi) + A 7 ($t$ 2 )]($tc i , 2 ) + h.c.| , (7) 

with the complex Higgs-doublet fields acquiring the vacuum expectation values v\ and V2- After 
diagonalising the mass matrix one obtains the physical Higgs fields h, H, A, and the Goldstone 
bosons G , G. We are interested in the Higgs self couplings, which can be written in terms of the 
dimcnsionless parameters Aj, i — 1, . . . , 7, appearing in ([7]). For the quartic couplings, one has 

^hhhh = 3 cos 4 aA 2 -I- 3 sin 4 aAi + 6 cos 2 a sin 2 a (A3 + A4 + A5) 

— 12 cos 3 a sin aA 7 — 12 sin 3 a cos a\§ , 

3 

^Hhhh = — 3 cosa sin 3 aAi + 3 cos 3 a sinaA 2 — - cos 2a sin 2a (A3 + A4 + A5) 
+3(3 cos 2 a sin 2 a — sin ci)\q + 3(cos 4 a — 3 cos 2 a sin 2 a)A7 , 

to be multiplied by Aff/v 2 . In the MSSM, the Higgs sector is constrained such that only two of the 
seven input parameters are free. Choosing tan/3 and Ma = Xqv 2 as basic input, one recovers the 
MSSM values 

Khhh = 3 cos 2 2a , 

^Hhhh — 3 cos 2a sin 2a . (8) 

As is well known, radiative corrections in the Higgs sector are large. The largest ubiquitous correction 
is given by 3G p-m 4 / (V2tt 2 sin 2 (3) ■ ln(m 2 /m 2 ). In our numerical analysis for the MSSM all one-loop 
and leading two-loop corrections to Higgs masses and couplings are included [36,37]. For a discussion 
about reconstructing the Higgs potential in the SUSY case, see [38]. 

3 Calculation 

From a computational point of view loop amplitudes with five or more external legs are challenging 
due to their combinatorial complexity. We only sketch here our method, which leads to a numerically 
stable algebraic representation of the amplitude. More details can be found in [39-44]. 

3.1 Structure of the amplitude 

The production processes pp — > hh and pp — > hhh proceed at the parton level via gluon fusion in 
combination with a quark loop. Here h stands for the Higgs boson of the SM or the light CP-even 
Higgs boson of the MSSM. Squarks in the loop have been neglected. Their contributions vanish for 
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large squark masses, due to their decoupling property. In the SM, only the top quark loop leads to 
a non- negligible cross section at the LHC. For triple Higgs boson production, 

g(pi, Ai) + g(p 2 , A 2 ) -> ft(p3) + MP4) + h(p 5 ) , (9) 

with ^ defining the 4-momenta and A1.2 specifying the gluon helicities, the Feynman diagrams are 
classified in Fig. [T] On the diagrammatic level 3-, 4- and 5-point topologies can be distinguished. 
Each topology involves a different combination of coupling constants. Here, we only discuss the 




Figure 1: Illustration of different topologies to the process gg — ► hhh. 

structure of the SM amplitude. In 2HDMs one has additional tree structures from the heavy Higgs 
boson H attached to the quark loop. The pentagon topology P contains no Higgs self coupling. The 
box topologies B are proportional to A3 , whereas the triangle topologies T\ and T2 are proportional 
to A4 and A§, respectively. Subsequently, the amplitude can be expressed as 

T(gg -» hhh) = S ab T R ^ e liM e v MT , (10) 



MT = X\ th M% + A 3 X 2 tth M% + A 4 X tth M^ A + Xl Xuh M% 2 . (II) 

The scattering tensor M pv can be decomposed in terms of metric tensors and external momenta 
and, using momentum conservation, be expressed as 

= AgT+ E ^iPjPl ■ (12) 
3,1=1,4 

By solving the Ward identities M £lP2 = 0, M PlE2 = 0, or equivalently using axial gauge conditions 
(e.g. £\ ■ P2, £2 'Pi), one can achieve a manifestly gauge invariant representation of the amplitude. 
Introducing the Abelian part of the gluon field strength tensor T^ v = e^Pj — Pj it reads 

M £ie2 = A tr(^^ 2 ) + £ Bji p2-Ti- Pj PX-T2-PL (13) 

J,J=3,4 
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The amplitude coefficients A, Bji arc equal to A, Bji, up to trivial factors. Bosc symmetry of 
the gluons and Higgs bosons leads to additional relations between them. After determining all 
amplitude coefficients, verifying the Ward identities and Bose symmetry served as a powerful check 
of our calculation. 

It is useful to decompose the amplitude further into helicity components. Due to parity invariance 
only two helicity amplitudes have to be known: 

M = M ++ , 

M~ + = M + ~ . (14) 
Applying spinor helicity methods [45,46], the polarisation vectors for ± helicities are given by 

[21] tr"(lz/2/i) 



H - v _ 



(12) 2s 12 ' 

(2-H1-) (2>|1- 
V2(21) V2[12] 



(15) 



with tr (lu . . . ) = [ti^ff" . . . ) — tr(75^/iy !y . . . )]/2 and the spinor inner products (ij) = (p i \pj), [ij] = 
(pt\pj), where \pf) is the Weyl spinor for a massless particle with momentum pj. 
This implies 



(12) V 2si2 2s 
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tr-(1324) tr"(1424) 



2S12 2S12 



43 ~ -B44 



i) (16) 



and 



(1 1 3 j 2 ) V 2si2 2si2 / 

<2— 14|1— > ^ tr-(1423) D t tr"(1424) 



(1-|4|2-) V 2si2 



B« + H3^ B44 ). (17) 



Contrary to the ++ case, it is not possible to factor out a global spinorial phase in the H — case 
without introducing denominators that in general aggravate numerical problems. 



3.2 Evaluation of the amplitude coefficients 

Our goal was the analytical reduction of all diagrams, to allow for algebraic cancellations of numer- 
ically dangerous denominators in the amplitude. These denominators are so-called Gram determi- 
nants which are induced by reduction algorithms of Lorentz tensor integrals. After generating all 
diagrams, using the QGRAF [47] program, we used FORM 3.1 [48] to perform the gamma matrix 
algebra and to project the diagrams on the helicity components and amplitude coefficients. Fur- 
ther, by applying the reduction algorithms for scalar and tensor integrals described in [42,43], we 
expressed all amplitude coefficients as a linear combination of scalar integrals. As scalar integral 
basis we chose 2-, 3- and 4-point functions (s^ = (pi +pj) 2 ): 

lt n {s l3 ,m 2 q ,m 2 q ) , 

li =i {sij , Skh Spr,m 2 q ,m 2 q ,m 2 q ) , 

lf =6 (sij , s k i, s pr ,ml,ml,m 2 q ,m 2 q ) , 

which were evaluated using LoopTools-2.2 [49]. The spurious UV pole of the 2-point integral cancels 
when adding all diagrams. The full amplitude is composed out of 12, 24 and 31 different 2-, 3-, 
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and 4-point functions. The complexity of the expressions is induced by the number of independent 
scales, which is seven here. One may chose S12, S23, S34, s 45 , S15, rv? q , m\. The coefficient of each 
function was exported to MAPLE to apply simplification algorithms. Schematically 

M x ^ =^simphfy[C fc AlA2 ]/ fe , I k e {/ 2 n , I 3 4 , 1 6 4 } . (18) 
k 

In this way, we could achieve expressions with a simple denominator structure allowing for a stable 
numerical evaluation. In the equal helicity case, Ai = A2, all Gram determinants cancel. In the 
opposite helicity case, Ai = — A2, one Gram determinant survives. The simplified expressions were 
then exported to Fortran code. Each of these steps was completely automatized. 

3.3 Numerical implementation 

In order to compute numerical results for hadron colliders, the differential partonic cross section has 
to be convoluted with parton distribution functions (PDFs) and integrated over the 2 — > 3-particlc 
phase space. We employed the gluon density of the MRST2002nlo PDF set [50], as implemented in 
LHAPDF [51], which also provides the strong coupling constant as function of the renormalisation 
scale. In the MSSM case, where the heavy CP-even Higgs boson H can be resonant, we used multi- 
channel Monte Carlo (MC) integration techniques [52,53] with phase space mappings based on [54,55] 
and the adaptive MC integration package BASES [56]. The relevant quartic Higgs couplings were 
implemented in the program HDECAY [37], which incorporates the routine FeynHiggsFast [57], in 
order to evaluate the radiatively corrected Yukawa- and Higgs couplings and also the Higgs widths, 
as discussed above. 



4 Results 

In this section we present and discuss the LHC cross sections for 2- and 3-Higgs boson production. 
We use the following parameters throughout: 

a s (M z ) = 0.120 
a(0) = 1/137.036 
m t = 178 GeV 
m b = 4.7 GeV 
m w = 80.41 GeV 
m z = 91.1875 GeV 

(19) 

For 2-Higgs (3-Higgs) boson production, the factorisation scale \xf and renormalisation scale [1r 
were set to fip = Hr = 2m h (3rajf). The strong coupling constant a s was taken at /Ur, but for the 
fine structure constant we used a(0). All results have been calculated with a MC error of 0.5% or 
less. 

In the following subsections we present SM cross sections for 2- and 3-Higgs boson production 
and describe how BSM scenarios allow for observable enhancements of the SM rates. 



4.1 Multi-Higgs boson production in the SM 

We begin with the gluon-fusion production cross section for 2 and 3 Higgs bosons. Note that in both 
cases large next-to-leading order corrections are expected, leading to iVfactors as large as 2 like in 
the case of single Higgs boson production [58,59], since the infrared structure of these processes is 
identical (with a large contribution from soft gluon effects). 
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m H [GeV] 

Figure 2: Total Higgs pair production cross section vs. itih at the LHC, as well as the individual 

cross sections for box (B) and triangle (T) topologies with equal (+ + / ) and opposite (+ — / — +) 

helicity components. 

In Fig. [21 we display the Higgs pair production cross section vs. run at y/s — 14 TeV. Our results 
show good agreement with [23] when PDF, scale and parameter uncertainties are taken into account. 
The total cross section falls from about 30 to 6 fb in the Higgs mass range from 100 to 200 GeV. In 

addition to the total cross section, the equal + + / and opposite H — / — h helicity components 

of the cross section are also shown. The triangle topologies only allow for a L = S = interaction, 
i.e. the H — / — h helicity component is zero. Overall, the opposite helicity component is more than 
an order of magnitude suppressed. Furthermore, a destructive interference effect is visible between 
the box and triangle topologies. 

In Fig. [3l the total cross section for triple Higgs boson production at the LHC is plotted vs. m#. 
Our results agree with the recent calculation of [26] within MC errors. Again one finds that the 
opposite helicity components of the cross section are more than an order of magnitude suppressed. 
In the figure, the pentagon, box, and triangle contributions are also shown separately. The latter, 
being proportional to A3 and A4, are suppressed relative to the box and pentagon topologies. Due to 
interference effects the contribution of the quartic Higgs self coupling to the total cross section is not 
negligible: it varies between +1% (tuh = 100 GeV) and -57.5% (m# = 200 GeV). The destructive 
interference pattern between triangle, box and respectively box and pentagon contributions is well- 
known. It can be understood from from the fact that the effective two-gluon n-Higgs boson operators 
contain a factor (—1)" [22]. As the self couplings increase with increasing Higgs mass, the box and 
triangle topologies become more and more important relative to the pentagon contribution. 

In Table Q] we give predictions for different values of the Higgs mass for the LHC and a 200 TeV 
Very Large Hadron Collider (VLHC). We also note that a change in the top mass best fit value from 
rrif = 178 GeV to 172.5 GeV leads to a 15% decrease of the cross section. 

In Table El we compare the relative importance of the different topologies and helicities at the 
VLHC. 

Furthermore, the gluon fusion cross section is proportional to a 2 s and thus very sensitive to 
renormalisation scale variations, as exemplified in Fig. 3] for fi = fj,R = fip. Here, the scale /x is 
varied around the central choice fio — 2 m# and fiQ — 3 m# for 2- and 3-Higgs boson production, 
respectively, by a factor /i//io G [1/3,3]. From this we estimate a scale uncertainty of about 50%. 
Thus, large if-factors due to higher order effects can be expected. In Table [3] we study the scale 
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Figure 3: The total 3-Higgs boson production cross section vs. ttih at the LHC. Equal and opposite 
helicity components of contributions from various topologies are shown. 





m H 


[GeV] 


120 


150 


180 


LHC 


a 


[fb] 


0.0623 


0.0267 


0.0126 


VLHC 


a 


[fb] 


9.55 


4.89 


2.98 



Table 1: Typical cross sections for triple Higgs boson production at the LHC and a 200 TeV VLHC. 



dependence of the cross section for 3-Higgs boson production by varying the renormalisation and 
factorisation scales independently for tuh — 120 GeV. Note that the gluon luminosity decreases with 
increasing scale /j,f, because in 3-Higgs boson production the momentum of the gluons has to be 
relatively high. Hence, the cross section shrinks with increasing [ip and The table demonstrates 
that varying fiR and hf in the same direction yields a conservative estimate of the scale uncertainty. 

Many years ago multi-Higgs boson production via gluon-fusion was studied in [60] in the heavy 
top limit. For single Higgs boson production via gluon fusion this limit is well known to be a good 
approximation [61]. In the context of multi-Higgs boson production, the heavy top limit has been 
applied at the leading [60] and next-to- leading level [62]. In [16], the quality of the heavy top 
approximation has been studied for Higgs pair production, and agreement at the O(10%)-level for 
the total cross section, but large discrepancies for kinematic distributions have been observed when 
comparing results for mt — * oo and physical m t . In Fig. [5] we compare our physical-m t results with 
results in the heavy top limit. While there is reasonable agreement in the 2-Higgs case for small 
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m H = 120 GeV 
m H = 200 GeV 


9.55 
1.93 


21.82 
3.97 


10.04 
4.76 


0.111 
0.129 


0.189 
0.262 


0.589 
0.125 


0.169 
0.163 



Table 2: Contributions of different topologies to triple Higgs boson production at a 200 TeV VLHC. 



8 




o-[10-*fb] 


fi F = m H 


3itih 


9lTlH 


fJ-R = m H 


9.71 


8.61 


7.66 


3iriH 


7.21 


6.38 


5.68 


9m,H 


5.57 


4.93 


4.39 



Table 3: Renormalisation and factorisation scale dependence of the gg — > HHH cross section for 
m H = 120 GeV at the LHC. 



Higgs masses, the heavy top limit fails completely in the 3-Higgs case. To better understand this 
observation, we study the variation of the cross sections with the internal quark mass m q for a fixed 
value of mH — 120 GeV in Fig.[6j We see that in the 2-Higgs case the heavy top limit accidentally 
agrees with the result for m q = mt (indicated by the vertical line). However, asymptotically the 
result for finite m q approaches the m q — > oo limit only for masses around 3 TeV. The same holds for 
the 3-Higgs case, but here the result for m q = mt is an order of magnitude smaller than the heavy 
top limit. In both cases the dominant contribution to the cross section comes from the kinematic 
regime close to the top pair threshold Sy ~ 4m \. We conclude that the heavy top limit is not 
applicable when calculating multi-Higgs boson production cross sections and should not be used in 
experimental studies. 

4.2 Multi-Higgs boson production beyond the SM 

If one allows for higher dimensional operators in the Higgs sector the trilinear and quartic Higgs 
self couplings are no longer directly related to the Higgs mass. This has motivated us to study the 
variation of the cross section with respect to A3, A4 and ran- In Fig. [7] we illustrate the variation 
of the cross section for 3-Higgs boson production with A3 and A4 {run is fixed), and in Fig. [8] the 
variation with A3 and to# (A4 is fixed). 

The variation of the cross section with A3, A4 and run is mainly due to phase space and interfer- 
ence effects. When m# increases, the phase space is reduced and the PDF and a s are to be taken 
at a higher value of x and a larger scale. All effects conspire and lead to a smaller cross section. 
The dependence on A3 is mainly due to the interference pattern. In Fig. [9] a slice of the contour plot 
Fig.[8]is shown for run = 160 GeV. When the trilinear coupling A3 becomes negative, the destructive 
interference between box and pentagon topologies turns into a constructive one, which explains the 
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Figure 6: The total cross section for 2-Higgs (left) and 3-Higgs (right) production vs. m q at the 
LHC. The m q — > oo limit is shown as horizontal line. The vertical line indicates the actual value of 
the top quark mass. 



rise of the cross section in this regime. Increasing the coupling to positive values beyond the SM 
value first enhances the destructive interference effect, but finally the box topologies win and the 
full cross section grows again. The minimum moves to lower values of A3 when the Higgs mass gets 
larger, due to the growth of the Higgs self-couplings. We note that the same effect happens in the 
2-Higgs case. It implies that even a precise measurement of the cross section alone would not lead 
to a unique determination of the trilinear coupling. 

This reasoning allows now to understand the variation of the cross section with A4 in Fig. [7J For 
negative values of A3 the box and pentagon topologies interfere constructively and the relative im- 
portance of the triangle contribution proportional to A4, A3 is further reduced, resulting in an almost 
fiat dependence on A4. For positive A3 the contribution of the triangle topologies is pronounced by 
destructive interferences between the various topologies leading to a slight variation with A4. As can 
be seen from Fig. [3 for A3/A§ M in the range 0.5 to 1.5 and A4/A4 in the range -3 to 3 the cross 
section varies from 0.03 to 0.1 fb. 

In principle the couplings are restricted only by unitarity bounds. To illustrate the effects of 
large couplings that approach the non-perturbative regime, we list in Table [5] cross sections for 
values I A3/ v| , I A4 1 ~ 47r. We see that non-perturbative effects in the Higgs sector may well lead to 
sizable triple Higgs cross sections of up to 30 fb at the LHC. 
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Figure 7: Contour plot showing the variation of the cross section for 3-Higgs boson production with 
A3 and A4 for ma = 120 GeV at the LHC. The numbers denote the cross section normalised to 
o-smOh = 120 GeV) = 0.0623 fb. 



As already pointed out in Section 2, in 2HDMs one finds two amplification effects for multi-Higgs 
boson production. Both effects are illustrated in Fig. [Til] for the MSSM. For small tan (3 = 3 the 
heavy CP even Higgs boson may become resonant in some of the topologies sketched in Fig. [JJ This 
is illustrated in Fig. [TOa . With the given choice of parameters one sees that for m/j > 109 GeV the 
H — ► hh channel opens up inducing a resonant amplification of the box and triangle topologies, B 
and T2. For rrih > 111 GeV also the H — > hhh channel opens up leading to an enhancement of the 
triangle topology T\ proportional to the quartic coupling Xnhhh- Since this triangle contribution is 
suppressed relative to the other topologies the effect is hardly visible in the plot (tiny peak next to 
the right vertical line). We see that in BSM scenarios resonant amplification may lead to triple Higgs 
production reaching a few fb (not taking into account the expected if- factor of 2), which would be 
observable at the SuperLHC or even at the LHC. 

For large tan/3 = 50 the top contribution is largely suppressed and the bottom loops become 
dominant due to the enhanced Yukawa couplings. Resonance effects on the other hand are now 

40 
32 
24 
16 
8 
4 
2 
1 

0.5 
0.1 



m H [GeV] 

Figure 8: Contour plot showing the variation of the cross section for 3-Higgs boson production with 
A 3 and m H for A 4 = Af M at the LHC, normalised to <7 S M{m H = 120 GeV) = 0.0623 fb. 
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Figure 9: Variation of the cross section for 3-Higgs boson production with A3 for uih — 160 GeV 
and fixed A 4 = A| M at the LHC, normalised to a S M(m H = 120 GeV) = 0.0623 fb. 



CT [fb] 


A 4 = -47T 





47T 


X3/V = — 4-7T 


28.0 


30.7 


33.4 





0.169 


0.0271 


0.0428 


47T 


12.2 


14.0 


15.8 



Table 4: 3-Higgs boson production cross sections for extreme choices of the Higgs couplings for 
m H = 150 GeV at the LHC. 

negligible, because the Higgs width Th is a factor 50 larger than at small tan j3 = 3. Based on the 
double Higgs production results in [21], we expect therefore that the LO triple Higgs production 
cross section in bottom quark fusion is smaller than in gluon fusion. As the pentagon contribution 
to the cross section is proportional to A?,t, all other topologies are much less relevant. For uia — » 00, 
m/j approaches its maximal value. In this limit one obtains the SM value for the cross section as 
shown in Fig. 110b . The minimum around ~ 115 GeV arises due to the interplay of the Yukawa- 
and Higgs self coupling. Although the cross section rises for decreasing m,h, it is always well below 
1 fb in the range considered. Note that in the 2-Higgs boson production case the bottom loops are 
more pronounced than in the 3-Higgs case, because the scalar loop integrals are probed in different 
kinematical regions. 

In principle, also squark loops have to be considered to obtain a complete prediction for hhh- 
production in the MSSM. In the heavy squark limit this contribution decouples and approaches zero 
in contrast to the quark contribution. For the present illustration of amplification effects we have 
neglected squark effects. 
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Figure 10: Cross section for 3-Higgs boson production vs. run at the LHC in the MSSM for tan j3 = 3 
(left) and tan/3 = 50 (right) including mixing effects (A t = ITeV, ^/TeV = -1(+1) for tan/3 = 
3 (50)). The vertical lines in the left plot indicate the 2- and 3-Higgs boson thresholds (H — ► hh and 
H hhh). 
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5 Summary 



We presented our calculation of the loop-induced processes gg — > HH and gg — > HHH, and 
discussed the resulting cross sections and their experimental accessibility at the LHC in the SM and 
beyond. The contributions from pentagon, box and triangle topologies exhibit strong interference 
patterns and large differences in the equal and opposite gluon helicity components. Generally, the 
opposite helicity component, which corresponds to a gluon pair with helicity 2, is suppressed by 
more than an order of magnitude. Furthermore, the triangle topologies are suppressed relative to 
the box and pentagon topologies. This results in a complex dependence of the cross sections on 
the trilincar and quartic Higgs self couplings. The cross section for 3-Higgs boson production varies 
strongly with the quartic Higgs self coupling. For to# = 100 GeV and 200 GeV it affects the cross 
section by +1% and —57.5%, respectively. This effect is due to the fact that the Higgs self coupling 
contributions are enhanced for higher Higgs masses and that there is a strong destructive interference 
between box and pentagon contributions. The estimated renormalisation and factorisation scale 
uncertainty of about 50% (variation by a factor 3) is typical for leading-order QCD cross sections. 
We also investigated the applicability of the heavy top quark approximation for 2- and 3-Higgs boson 
production and demonstrated that the heavy top limit is not applicable when calculating multi-Higgs 
boson production cross sections. In summary, we find that the SM cross sections for 3-Higgs boson 
production are too small for observation at the LHC. The measurement of the SM trilinear Higgs 
self coupling in 2-Higgs boson production, on the other hand, may be feasible at a luminosity 
upgraded LHC, termed SuperLHC, collecting 6000 fb -1 of data, as long as QCD backgrounds are 
not prohibitively large [8,63]. 

The experimental prospects improve if favourable extensions of the SM are realized in nature. We 
have demonstrated that cross sections for triple Higgs boson production can be as large as O(10) fb 
if one allows for higher dimensional operators or considers 2-Higgs-doublet models, e.g. the MSSM. 
Here, two amplification effects have been analysed. First, we studied Yukawa coupling enhancements 
through mixing that lead to important bottom loop contributions. For the 3-Higgs case the cross 
section remains below 1 fb for the parameters considered and tan (3 = 50. Secondly, we demonstrated 
that resonance effects due to an internal heavy Higgs boson decaying into 2 or 3 light Higgs bosons can 
enhance the 3-Higgs cross section into the potentially observable O(fb) region. Higher dimensional 
operators lead to essentially unconstrained trilinear and quartic couplings. Strong coupling of the 
order | A3 /u | , | A4 1 ~ A-k leads to cross sections of up to 30 fb. This implies that already the LHC 
should be able to restrict the A3-A4 plane. Although these bounds will not be very restrictive, 
they may still exclude parameter regions that exhibit non-perturbative effects, which would be an 
important qualitative finding. A more stringent bound on A3 alone will be obtained from Higgs pair 
production. We conclude that multi-Higgs boson production at the luminosity-upgraded LHC is an 
interesting probe of Higgs sectors beyond the SM and warrants further study. A 200 TeV VLHC 
would of course further improve the sensitivity. 
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